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Nematic droplets suspended in the isotropic phase of the same substance were subjected to alter- 
nating electrical fields of varying frequency. To keep the system at a constant nematic/isotropic vol- 
ume ratio with constant droplet size we carefully kept the temperature in the isotropic/nematic co- 
existence region, which was broadened by adding small amounts of a non-mesogenic liquid. Whereas 
the nematic droplets remained spherical at low (in the order of 10 Hz) and high frequencies (in the 
order of 1 kHz), at intermediate frequencies we observed a marked flattening of the droplets in 
the plane perpendicular to the applied field. Droplet deformation occurred both in liquid crystals 
(LCs) with positive and negative dielectric anisotropy. The experimental data can be quantitatively 
modelled with a combination of the leaky dielectric model and screening of the applied electric field 
due to finite conductivity. 
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I. INTRODUCTION 

Liquid crystal droplets dispersed in isotropic me- 
dia have attracted strong interest. Driven by new 
applications, during the last few years the research 
focus has shifted from studying the director field 
and defect structure in quasi-static situations (see, 
e.g.^ ^ El HI m [6|) to dynamic processes under 
the influence of external fields [3 [HI IS HO]. Due 
to the anisotropic nature of the liquid crystalline 
phases, these droplets show a much richer physics 
than purely isotropic systems. 

The literature on the behavior of liquid crystalline 
droplets subjected to external fields is dominated by 
studies on droplets suspended in a non-mesogenic 
carrier fluid. Owing to the increasing importance 
of polymer dispersed liquid crystals in various ap- 
plications, the scientific interest in the response of 
droplets to external field has grown significantly (for 
reviews see, e.g., [H H [El IIIIIS] ) • In addition, 
magnetic [2 and shear fields [9^ have been applied 
to these systems. Only few studies have been car- 
ried out on liquid crystals in the coexistence region. 
Park and coworkers [7 studied nematic droplets in 
a quasi 2-dimensional system submitted to an elec- 
tric field superimposed with a temperature gradi- 
ent. Dolganov and coworkers [6] investigated ne- 
matic droplets confined by a surrounding smectic 
phase. 

Inevitable impurities in liquid crystal samples are 
known to open up a coexistence region between the 
nematic and isotropic phase (see, e.g., [161 [13 UHl 
[T9l [20]). This coexistence region allows to have 
nematic droplets immersed in the isotropic phase 



of the same liquid at constant sample temperature. 
The physics of droplets in the isotropic-nematic co- 
existence region is peculiar because, unlike typical 
liquid-liquid interfaces, the interfacial tension is very 
low. Since droplet deformation always implies an 
increase in surface area, the higher the surface ten- 
sion, the more rigid the droplets appear. In contrast, 
droplets in the nematic-isotropic coexistence region 
have a surface tension [211 [221 [23 [21] of the order 
of 10~^ Nm and can be deformed very easily. 

In this paper, we discuss electric field-induced 
droplet deformations that are peculiar to systems 
with low interfacial tension. We put special empha- 
sis on the influence of field frequency, and the actual 
electric field strength. Moreover, we took special 
care to avoid electrochemical reactions in the sam- 
ple that might influence the results. The paper is 
organized as follows: In Sec.[n]we describe the ma- 
terials and the experimental set-up. In Sec. [Ill A} we 
present the results on the droplet deformation in the 
colloid- free systems. Sec. |IIIB| is devoted to inves- 
tigating the underlying hydrodynamic flow field by 
adding tracer colloids to the system. We then com- 



pare our results to theoretical predictions in Sec. IV 



II. EXPERIMENTAL DETAILS 

Both the electric and dielectric properties of 
nematic liquid crystals are markedly anisotropic. 
Whereas the conductivity a typically exhibits posi- 
tive anisotropy 



crii > cr^ 



(1) 



^Electronic address: auhammer@mpip-mainz.mpg.de| 



(the indices || and ± indicating the conductivity par- 
allel and perpendicular to the nematic director n, 
respectively), the anisotropy of the dielectric con- 
stant e can vary (and change sign) between different 
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FIG. 1: Schematic phase diagram of a Uquid crystal with 
added impurities (in this work: octane). With increas- 
ing octane concentration, the nematic isotropic phase 
transition shifts to lower temperatures and a nematic- 
isotropic coexistence region opens. The temperature Te 
and the octane concentration Ce (gray lines) were chosen 
be well within the coexistence region. 
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FIG. 2: The sample was placed between two ITO coated 
glass slides (purchased form SPI Supplies) without di- 
rect contact between the sample and the electrodes. The 
sample thickness is approximately 180 fim and the elec- 
trode spacing 500 /xm. After loading, the sample was 
sealed with UV glue (Norland Products, USA). Within 
the experimental time, typically less than 5 hours, no 
changes in the properties of the sample were observed. 



substances 

(2) 

We chose a typical representative for either 
class: methyloxybenzylidenbutylaniline (MBBA, 
Frinton Laboratories, USA, Ae < 0) and pentyl- 
cyanobiphenyl (5CB, Chemos, Germany, Ae > 0). 
Both chemicals were used as supplied. Impuri- 
ties like polymers or low molecular weight organic 
sovents in liquid crystals lead to coexistence of the 
nematic and the isotropic phase within a certain 
temperature interval [ifl HH |25l |26] (see Fig{l] for 
a schematic representation of the phase diagram). 
In our experiments, adding octane (Aldrich, ana- 
lytic grade) at a concentration of Cg = 2 vol% re- 
sulted in a coexistence interval of the order of 1.5 K, 
which is sufficient to override the effect of any im- 
purities inevitably present in commercially available 
liquid crystals. The results presented in this paper 
did not change for octane concentrations between 
1 and 4 vol%. Poly (methyl methacrylate) particles 
(PMMA) with a radius of approximately 500 nm 
sterically stabilized by chemically grafted poly (12- 
hydroxy-stearic acid) molecules [27l [28] were used 
as tracer particles. 

As depicted in Fig. |2j the sample was placed be- 
tween two glass slides of 170 /im each. The outer 
side of these slides was covered with transparent 
indium-tin-oxide (ITO) electrodes. This configura- 
tion ensured that there was no electric contact be- 
tween the electrodes and the sample. The surface of 
the glass slides was otherwise untreated. The sam- 
ple cell was temperature controlled to a precision 
better than 0.05 K using a Linkam heating stage 
(THMS 600). All observations were performed using 
an Olympus BX51 optical microscope with a long 
working distance 20x objective and equipped with a 



CCD camera. For better contrast and good resolu- 
tion of droplet edges, the sample was placed between 
parallel polarizers. The ITO electrodes allowed to 
subject the sample to an alternating electric field 
parallel to the light path of up to 0.4 V//im and 
with frequencies up to 5 kHz. An oscilloscope was 
used to ensure that the applied voltage remained 
constant over the whole frequency range. 

After placing the sample in the measurement cell 
we first thermally equilibrated the system for 20 to 
30 minutes at a temperature several degrees above 
the isotropic/nematic phase transition temperature 
Tin. Then we slowly cooled the system into the co- 
existence region (typically 0.1 to 0.2 K below Tjn)- 
Until the diameter of the nematic droplets formed 
in the sample reached 50 to 150 jam. The sample 
was then equilibrated at this temperature for 10 to 
15 minutes. In this way, the temperature could at- 
tain a constant value throughout the sample and 
the octane concentrations in both phases could ad- 
just to the values given by the phase diagram. Since 
the density of nematic droplets is a few 10~^ g/cm^ 
greater than that of the isotropic matrix (see [29] 
and references therein), they sink to the bottom of 
the cell during this second relaxation time. 

Three types of frequency dependent measure- 
ments were performed: i) logarithmic increase of the 
frequency with time from low to high frequencies, 
i.e., constant time per frequency decade, ii) loga- 
rithmic decrease from high to low frequencies, and 
iii) switching directly from zero field to the desired 
frequency and field strength. In all three cases, the 
two-dimensional projection of the droplet size onto 
the observation plane was measured as a function of 
time using ImageProPlus software (Media Cyber- 
netics, USA). 
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FIG. 3: a) Images of nematic droplets of MBBA and 5CB dispersed in the isotropic phase in the absence of an 
appUed electric field, b) In the series of images taken during frequency sweeps (sweep time: 100 s for MBBA: 10 Hz 
to 1 kHz and for 5CB: 50 Hz to 5 kHz, voltage 0.4 V//im) the droplet shape and apparent size changes as a function 
of frequency. The scaling bar indicates 100 /xm. 



III. RESULTS 
A. Colloid-free system 

Figure |3] shows the results of typical frequency 
sweeps in MBBA and 5CB. Part a) shows field- free 
droplets, whereas in part b) all images were taken 
at 0.4 V//im. Now, the director orientation in the 
droplets can be assumed to be predominantly per- 
pendicular (MBBA) or parallel (5CB) to the optical 
path. In the absence of any applied field, there was 
no significant difference between MBBA and 5CB 
droplets. The observed texture of the droplets was 
consistent with spherically symmetric director ori- 
entation inside the droplets [30]. 

When an electric field was applied, the texture of 
the MBBA droplets changed, as shown in Fig. [3] 
This was caused by the negative dielectric contrast 
Ae of MBBA which leads to a predominantly in- 
plane director orientation. In the case of 5CB, the 
texture hardly change when turning on a low fre- 
quency field. This might be surprising since many of 
the material parameters (like the elastic constants) 
of the two samples are similar (see Tab.|l|. In 5CB, 
however, the applied field did not change the sym- 
metry along the optical path. Turning the director 
inside of the droplet parallel to the field and to the 
optical path keeps the rotational symmetry of the 
droplet projection through the microscope. This is 
different for MBBA, where the electric field turns 
the director perpendicular to the optical path and 
thus breaks the rotational symmetry of the projec- 
tion image of the droplet. 

Changing the frequency of the electric field not 
only changed the apparent size of the droplets, but 
also their internal texture, indicating a change of di- 
rector orientation. In 5CB, this effect increased with 
frequency, corresponding to a good director orienta- 
tion parallel to the field. In MBBA, the change was 



not visible due to the negative dielectric constant. 
From the evolution of the texture we can conclude 
that the higher the frequency,the better the orienta- 
tion inside the droplets, ie., the higher the electric 
field inside the sample. 

All commercial liquid crystals have a small but fi- 
nite conductivity based on ionic impurities. Since 
the electrodes are electrically insulated from the 
sample, we expect that the field in the sample is 
smaller than the applied electric field due to screen- 
ing by ion migration. At low frequencies, the ions 
can easily follow the applied field and screen the 
sample, whereas at sufficiently high frequency, the 
limited ion mobility prevents this. Thus, the field 
is more markedly reduced at low than at high fre- 
quencies. The increase of the field as a function of 
frequency is responsible for the increasing orienta- 
tion of the nematic droplets in 5CB as shown in 
Fig. Isb). In Sect. IV A we shall give a quantitative 



description of this effect. 

As can be seen in Fig. ^p) (5CB) and in Fig. [4] 
(MBBA), the amplitude of the diameter change de- 
pends on the initial droplet size. In both 5CB and 
MBBA, droplets with an initial (field-free) diame- 
ter of 50 - 200 /im exhibited marked changes of the 
apparent size after field application. In MBBA, the 
amplitude ratio, i.e., the ratio of the maximum di- 
ameter to the initial diameter, varied from 1.2 for 
droplets initially 50 /im in diameter to values up to 
2 for initial diameters of 140 /im and decreased again 
for even larger droplets. A similar behavior was ob- 
served in 5CB. Note that the frequency for the max- 
imum diameter in 5CB is higher by about one order 
magnitude than in MBBA (see also Sect. IV). 



Despite all these differences, the observed behav- 
ior was very similar in both systems. Furthermore, 
the samples could be used in many frequency sweeps 
(at least for several hours) without any noticeable 
change in the frequency response. The electrical in- 
sulation between the electrodes and the sample sup- 



3 



G. K. Auernhammer, et al. 



Quantity 


MBBA 

Value Ref. 


5CB 

Value Ref. 


Interfacial tension 7 


1.5 • 10"^ Nm"^ 23 


1.6 • 10"^ Nm"^ [22 


Dielectric constant e\\ 

e± 

^iso 

Dielectric ratio S = ^ 


4.85 31_ 
5.1 [31] 
5 [33] 
0.92 


15.5 [32 
7.5 [32] 
10.5 [32] 
0.65 


Electric conductivity am 
Conductivity ratio R = ^ 


5 • 10-^ Sm-i 

0.9 31_ 


1.5 • lO-'^ Sm-i 

1.1 [34 


Viscosity ratio m = ^ 


2.1 7 


0.7 [7 


Elastic constant Ki 
K2 
Ks 


6 • 10"^^ N 35 
3.8 • 10"^^ N 35 
7.5 • 10"^^ N 35 


6.2 • 10"^^ N [36 
3.9 • 10"^^ N [36 
8.2 • 10"^^ N [36 



TABLE I: Material parameters of MBBA and 5CB in the coexistence region. The electric conductivities were fitted 
to match the model with the experimental data as shown in Fig. [9] Still these values come within the limits typically 
expected for these materials. Note that S is calculated for oriented droplets in which the director is perpendicular 
(MBBA) or parallel (5CB) to the applied field, with Cm and Cd being the dielectric constants of the continuous 
medium and the droplets, respectively. 
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FIG. 4: The typical evolution of MBBA droplets while 
changing the frequency in three different ways (up: log- 
arithmic sweep from 10 Hz to 1 kHz; down: logarithmic 
sweep from 1 kHz to 10 Hz; jump: each frequency is 
directly addressed from a field free state) . The apparent 
droplet diameters are normalized by their diameter be- 
fore applying the field, which were measured to be 135 
/xm (up), 114 /im (down), and 91 /xm (jump). The dif- 
ferent amplitudes of the curves are due to the different 
initial droplet size. The sweeps in the major plot were 
performed with 100 s sweep time. The inset shows the 
time dependence of the relative diameter after turning 
on the voltage (jump at 100 Hz). 



presses electrochemistry as a possible mechanism of 
degradation. However, the sample properties did 
change slightly in samples stored for more than 
24 hours. For that reason, we only consider data 
taken within the first 5 hours after sample prepa- 
ration. This virtually excludes field-induced dam- 
age as has been reported in experiments on electro- 
hydrodynamic convection (see, e.g., ^37]). 



The shape change was almost independent of fre- 
quency history as shown for MBBA in Fig. [4j The 
amplitude variations are in the range expected for 
droplets of different size. However, the frequency of 
maximum amplitude fmax was slightly less in the 
case of decreasing frequency sweeps. This indicates 
a small memory effect which might arise from ohmic 
heating in the ITO layers or from relaxation phe- 
nomena as shown in the inset to Fig. [4j After a 
jump in the applied voltage it takes 10 - 15 s until 
the droplet has attained its stable shape. 

As temperature is kept constant, also the nematic 
volume fraction and droplet volume can be consid- 
ered constant. Therefore, the observed increase in 
apparent diameter must be caused by flattening of 
the droplets, i.e., the droplets take an oblate (disk- 
like) shape with symmetry axis along the electric 
field (as sketched in Fig. [4|. This behaviour differs 
from that predicted by purely dielectric considera- 
tions [38l. 



B. Tracing the flow fleld 

To test whether droplet deformation is of hydro- 
dynamic origin, we added tracer colloids to the sam- 
ple mixture prior to filling the chamber. In the co- 
existence region, these colloids remain exclusively in 
the isotropic part of the sample. Thus they visualize 
the flow field around the droplets. 

These experiments revealed a strong flow within 
the isotropic phase around the nematic droplets. 
The tracer colloids moved outwards close to the 
droplet surface, turned upwards well inside the 
isotropic areas, returned above the droplet and then 
closed their convection role by moving downwards 
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FIG. 5: Sketch of the flow field around the nematic 
droplets as observed from tracer colloids. The flow field 
is axisymmetric around the center axis of the droplet 
(dotted line). 
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IV. THEORETICAL CONSIDERATIONS 

A. Screening 

To describe the screening of the electric field due 
to the finite conductivity of the sample, we use the 
model circuit sketched in the inset in Fig. jT] Here, 
we model the glass slides as pure capacitances, Cgi^ 
and the sample as a capacitance, Cg, in parallel 
with an resistance, Rg (assumed to be ohmic). The 
impedance (Zceii) of this model circuit is given by 
the sum of the impedances of the glass slides, Zgi, 
and the sample, Zs 



Zgl 

z. 



Kg 



2Z, 



Zsj 



(3) 

(4) 
(5) 



with i = u is the angular frequency of the 

applied field, the capacities given by C = {eocA/d)^ 
the ohmic resistance Rs = d/{Aa)^ the sample area 
by A, and the thickness d = dgi = dg. The voltage, 
[/g, actually applied to the sample given by 



FIG. 6: Characteristic flow velocities around a droplet 
(filled squares) and the corresponding apparent diame- 
ter of the droplet (open squares). To have a velocity 
measure which is independent from the actual droplet 
size, all flow velocities are measured at the perimeter 
of the maximum droplet size, i.e., for this droplet of a 
maximum diameter of 75 /xm at about 40 /xm from the 
droplet center. 



from above the droplet to its side at some distance 
from the interface. The resulting toroidal flow field 
is sketched in Fig. |5j The projected velocity of the 
tracer colloids was highest close to the surface of 
the droplet. This indicates that the flow is driven 



by the interface (see Sect IV for more details). Via 
viscous friction, the toroidal flow field leads to a 
net outward force on the droplet which leads to the 
transition from the spherical to the oblate droplet 
shape. 



Uo 



^cell 



2ujeoes — 2ia + ooe^egi ' 



(6) 



where cq = 8.85 • 10~^^F/m is the permittivity of 
vacuum, Cs and Cgi are the dielectric constants of 
the glass slides and the sample, and a is the con- 
ductivity of the sample. Taking as typical values for 
MBBA = 5, cr = 0.5 • 10-'^5'/m, Cgi = 5, and 
Uq = 2 OOF, we found that the voltage actually ap- 
plied to the sample, [/g, reaches the value expected 
from purely dielectric considerations, Jjdzeiect. ^^^y. 
for frequencies of a few hundred Hz (see also Fig. [7| 



u: 



dielect. 



Uo 



2e., 



Uo 
3 



lim Us 



(7) 



As expected, the sample is essentially field-free at 
low frequencies 



lim Us 



0. 



(8) 



The velocity maximum at intermediate frequen- 
cies coincides with maximum deformation (Figj6|. 
Within the accuracy of the data, the flow velocity 
vanishes at high and low applied frequencies. Flow 
fields around droplets also lead to long-range inter- 
action between the droplets [39] . In our low surface 
tension system, this interaction leads to deviations 
from the circular shape, as shown in Fig. [sJd). 



B. Leaky dielectric model 

The "leaky dielectric model" describes the behav- 
ior of a droplet immersed in an immiscible medium 
when subjected to an electric field. In this model, 
both the droplet and the continuous phase are as- 
sumed to be slightly conducting dielectrics. The 
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FIG. 7: Modeling the effects of a small but finite con- 
ductivity of the samples. For typical values for MBBA 
(cs = 5, CT = cTm — 0.5 • 10~^S'/m, Cgi — 5, and 
Uo = 200V) the actual field in the sample reaches only 
at a few hundred Hz its purely dielectric value. 




Here, we include a brief review of the leaky di- 
electric model (see the references cited above for de- 
tails). Any electric current j induced by an applied 
field has to fulfill the continuity equation V • j = 0. 
Additionally, the electric field is submitted to the 
usual boundary conditions at the interface between 
the droplet and the medium V • eeoE = and 
V X E = 0. The combination of these conditions 
typically leads to a finite electrical charge density 
on the droplet surface. As illustrated in Fig. [Sj the 
interaction of these surface charges with the applied 
electric field induces electrokinetic effects, which ul- 
timately lead to a fluid motion parallel to the droplet 
surface. Depending on the ratios of the conductivi- 
ties, viscosities, and dielectric constants of both me- 
dia, the flow above and below the droplet is directed 
outwards or inwards. Outward (inward) flow leads 
to an oblate (prolate) deformation of the droplet. 
The deformation consists of a stationary and a time- 
dependent part. The amplitude of the latter turned 
out to be below the detection limit of our set-up, i.e., 
is negligible in our system. In the notation adopted 
by Saville and coworkers [42] [43], the stationary de- 
formation of the droplet Ds = {di — (i2)/(<^i + <^2) 
(with di and d2 being the principal diameters of the 
droplet) can be expressed by 



FIG. 8: The combination of electric boundary conditions 
and continuity equation lead to surface charges on the 
droplet surface. These surface charges interact with the 
applied electric field and induce fluid motion tangential 
to the droplet surface. 

foundations of that model have been laid by Taylor 
[40] . Torza Saville [42l|43], and their coworkers. 
Recently this model has experienced more detailed 
experimental [39l [44] and numerical [45l [46] verifi- 
cations. 



Dstat = T^^oem (9) 

16 7 

with a representing the droplet radius in the absence 
of an electric field and E the actual field in the sam- 
ple. Since the interfacial tensions 7 enters in the de- 
nominator, it is clear that the extremely small value 
of 7 ~ 1.5 • 10~^ N/m enables large deformation of 
the droplets. The function ^ determines the type of 
the deformation. It is positive (negative) for prolate 
(oblate) droplets 



S^RjU + 14M) + 15^^(1 + M) + ^(19 + 16M) + R^Stu;{M + 1){S + 2) 
5(M + 1)[52(2 + R)^ + i?V2cj2(l + 25)2] • 



(10) 



Here, R = (Jd/cTm denotes the ratio of the conductiv- 
ities, S = Cm/^d the ratio of the dielectric constants, 
M = Vrn/l-^d the ratio of the viscosities, r = eeo/cr 
the characteristic time, uo the angular frequency of 
the applied electric field, and indices d and m indi- 
cate the properties of the droplet and the continuous 
medium, respectively. <l> is a monotonic function of 
uo which may change sign as a function of cj, but has 
no maximum or minimum at finite uo. To obtain 
a maximum in the deformation we have to include 



the frequency dependent electric field in the sample 
discussed above. 



C. Comparison with our data 

The solid line in Fig. [9] combines our estimate of 
the actual voltage applied to the sample, Es = Us/d 
from Eq.(|6|, with the frequency dependent defor- 
mation given by Eq. (10). In this comparison we 
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FIG. 9: The deformation Ds = {di - d2) / {di + d2) (with 
di and d2 being the principal diameters of the droplet) 
as calculated from the data of Fig. |4] (assuming a con- 
stant droplet volume) is compared to the combination 
of screening and the leaky dielectric model (solid line). 
The modelling curve has been calculated with the ma- 
terial parameters given in Tab. [l] an applied ac voltage 
with 200 V amplitude, a droplet radius of 60 /im, and 
an electrode distance of 500 /xm. 
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FIG. 10: Frequency dependent deformation of nematic 
droplets in a 5CB/ionic liquid mixture (about 0.01 vol% 
of l-butyl-3-methylimidazolium trifluoromethansulfon- 
imid in 5CB). 



the samples (5CB) with an ionic liquid (1-butyl- 
3-methylimidazolium trifluoromethansulfonimide) . 
This strongly suppressed all field-induced effects. 
At a concentration of a few vol% neither a change 
in texture nor a deformation was observed. Even 
at much lower charge concentrations, both effects 
were detectable but partially suppressed compared 
to pure samples (see Fig. 10). We attribute this 



strong reduction to the high concentration of free 
ions which leads to a significantly higher conductiv- 
ity and thus to better screening. 



V. CONCLUSION 

The extremely low interfacial tension at the 
nematic-isotropic interface allows for easy deforma- 
tion of nematic droplets in the isotropic-nematic co- 
existence interval. 

The reported frequency-dependent oblate (disc- 
like) deformation of nematic droplets induced by 
an applied electric ac field extents the behaviour 
expected in "leaky" dielectrics. A part from the 
anisotropy in the material parameters, the nematic 
degrees of freedom seem not to play a significant 
role in mechanisms leading to the droplet deforma- 
tion. Working with insulated electrodes suppressed 
electrochemical reactions at the electrodes but also 
made it necessary to include screening in the the- 
oretical description. As characteristic for leaky di- 
electrics, an interface-driven, frequency-dependent 
hydrodynamic flow field around the droplets was in- 
duced by the electric field. The effective field how- 
ever was reduced due to screening by field-induced 
charges at the interface, caused by the finite conduc- 
tivity of the sample. Accordingly, combining screen- 
ing with the " leaky dielectric" model accurately de- 
scribes the experimentally observed droplet defor- 
mation. 



used the material parameters (see Tab. |l| given in 
the literature and chose the value of the only free 
parameter in the model, the conductivity, to fit 
the measured deformations. The resulting values 
for the conductivity {(Jmbba = 5 • 10~^ S/m and 
(^bCB = 1.6 • 10~^ S/m) are in the range expected 
for liquid crystals. With the above assumptions, the 
combination of leaky dielectric model and screening 
due to the finite conductivity accurately describes 
the measured data. 

To verify the effect of free charges, we increased 
the concentration of free ions by mixing one of 



Acknowledgment 

It is pleasure for GKA to acknowledge fruitful 
discussions with Harald Pleiner and Patricia E. 
Cladis. We are grateful to Jochen S. Gutmann for 
synthesis of the ionic liquid and to Gabi Schafer 
and Norbert Hohn for their help with the data 
treatment. This work has been supported by the 
Deutsche Forschungsgemeinschaft through the SEE 
TR6 'Physics of Colloidal Dispersions in External 
Fields'. 



[1] G. P. Crawford and S. Zumer (Eds.), Liquid Crys- 
tals in Complex Geometries (Taylor and Francis, 



London, UK, 1996). 
[2] S. Candau, P. Le Roy, and F. Debeauvais, Mol. 



7 



G. K. Auernhammer, et al. 



Cryst. Liq. Cryst. 23, 283 (1973). 
[3] J. W. Doane, N. A. Vaz, B.-G. Wu, and S. Zumer, 

AppL Phys. Lett. 48, 268 (1986). 
[4] P. K. Chan, Liquid Crystals 26, 1777 (1999). 
[5] T. C. Lubensky, D. Pettey, N. Currier, and 

H. Stark, Phys. Rev. E 57, 610 (1998). 
[6] P. V. Dolganov, H. T. Nguyen, G. Joly, V. K. Dol- 

ganov, and P. Cluzeau, Europhys. Lett. 78, 66001 

(2007). 

[7] C. S. Park, N. A. Clark, and R. D. Noble, Phys. 

Rev. Lett. 72, 1838 (1994). 
[8] B. L Lev, V. G. Nazarenko, A. B. Nych, D. Schur, 
P. M. Yamamoto, and H. Yukuyama, Phys. Rev. E 
64, 021706 (2001). 
[9] Y. Wu, W. Yu, C. Zhou, and Y. Xu, Phys. Rev. E 
75, 041706 (2007). 
[10] R. El-Sadek, M. Roushdy, and J. Magda, Langmuir 

23, 7907 (2007). 
[11] D. Coates, J. Mater. Chem. 5, 2063 (1995). 
[12] L. Bouteiller and P. LeBarny, Liq. Cryst. 21, 157 
(1996). 

[13] T. J. Bunning, L. V. Natarajan, V. P. Tondiglia, 
and R. L. Sutherl, Ann. Rev. Mater. Sci. 30, 83 
(2000). 

[14] S. Tjong, Mater. Sci. Eng. Reports 41, 1 (2003). 
[15] Y. J. Jeon, Y. Bingzhu, J. T. Rhee, D. L. Cheung, 

and M. Jamil, Macromol. Theory Simul. 16, 643 

(2007). 

[16] H.-W. Chiu and T. Kyu, J. Chem. Phys. 103, 7471 
(1995). 

[17] C. Shen and T. Kyu, J. Chem. Phys. 102, 556 
(1995). 

[18] Z. Lin, H. Zhang, and Y. Yang, Phys. Rev. E 58, 
5867 (1998). 

[19] S. DasGupta, P. Chattopadhyay, and S. K. Roy, 

Phys. Rev. E 63, 041703 (2001). 
[20] P. K. Mukherjee, J. Chem. Phys. 116, 9531 (2002). 
[21] D. Langevin and M. A. Bouchiat, Mol. Cryst. Liq. 

Cryst. 22, 317 (1973). 
[22] S. Faetti and V. Palleschi, Phys. Rev. A 30, 3241 

(1984). 

[23] R. Williams, Mol. Cryst. Liq. Cryst. 35, 349 (1976). 

[24] H. Yokoyama, S. Kobayashi, and H. Kamei, Mol. 
Cryst. Liq. Cryst. 129, 109 (1985). 

[25] F. Roussel, J.-M. Buisine, U. Maschke, X. Coqueret, 
and F. Benmouna, Phys. Rev. E 62, 2310 (2000). 

[26] B. Ullrich, E. Ilska, N. Hohn, and D. Vollmer, Prog. 
Colloid. Polym. Sci. 113, 142 (2006). 

[27] L. Antl, J. W. Goodwin, R. D. Hill, R. H. Ottewill, 
S. M. Owens, S. Papworth, and J. A. Waters, Col- 
loids and Surfaces 17, 67 (1986). 

[28] G. Bosma, C. Pathmamanoharan, E. H. A. de Hoog, 
W. K. Kegel, A. van Blaaderen, and H. N. W. 
Lekkerkerker, Journal of Colloid and Interface Sci- 
ence 245, 292 (2002). 

[29] J. Deschamps, J. P. M. Trusler, and G. Jackson, J. 
Phys. Chem. B 112, 3918 (2008). 

[30] R. Ondris-Crawford, E. P. Boyko, B. G. Wagner, 
J. H. Erdmann, S. Zumer, and J. W. Doane, J. 
Appl. Phys. 69, 6380 (1991). 

[31] D. Diguet, F. Rondelez, and G. Durand, C. R. Acad. 
Sci. Paris B 271, 954 (1970). 

[32] S. Chandrasekhar, Liquid crystals (Cambridge Uni- 
versity Press, Cambridge, UK, 1977). 



[33] P. G. de Gennes and J. Prost, The Physics of Liquid 
Crystals (Clarendon Press, Oxford, 1993). 

[34] J. Jadzyn and P. Kedziora, Mol. Cryst. Liq. Cryst. 
145, 17 (1987). 

[35] H. Ivan, The Journal of Chemical Physics 57, 1400 
(1972). 

[36] L W. Stewart, The static and dynamic continuum 
theory of liquid crystals (Taylor and Francis, Lon- 
don, 2004). 

[37] H. Richter, A. Buka, and L Rehberg, Phys. Rev. E 

51, 5886 (1995). 
[38] G. Taylor, Proc. R. Soc. Lond. A 280, 383 (1964). 
[39] J. C. Baygents, N. J. Rivette, and H. A. Stone, J. 

Fluid Mech. 368, 359 (1998). 
[40] G. Taylor, Proc. R. Soc. Lond. A 291 (1966). 
[41] S. Torza, R. G. Cox, and S. G. Mason, Phil. Trans. 

R. Soc. Lond. A 2969, 295 (1971). 
[42] O. Vizika and D. A. Saville, 239, 1 (1992). 
[43] D. A. Saville, Annu. Rev. Fluid Mech. 29, 27 

(1997). 

[44] N. Bentenitis, S. Krause, and K. Benghanem, Lang- 
muir 21, 790 (2005). 

[45] A. Ferndez, G. Tryggvason, J. Che, and S. L. Cec- 
cio, Phys. Fluids 17, 093302 (2005). 

[46] J. Zhang and D. Y. Kwok, J. Comput. Phys.. 206, 
150 (2005). 



8 



